NYO-1480-7 


Courant  Institute 

of  Mathematical  Sciences 

MAGNETO-FLUID  DYNAMICS  DIVISION 


MP -41 


Nyo-i48o-7 


Asymptotic  Equivalence 
of  the  Navier-Stokes  and 
Nonlinear  Boltzmann  Equations 

HAROLD  GRAD 

September  30,  1964 

AEG  Research  and  Development  Report 


NEW    YORK    UNIVERSITY 


\ 


Mzm 


»::,:'*'  Jl 


■r:i^r-^- 


j^'^'m,''^  ;^  '=**• 


•-  *-,  i>,.  »^^  *  ■'  B    « 


This  report  was  prepared  as  an  account  of  Government  sponsored  work.  Neither 
the  United  States,  nor  the  Commission,  nor  any  person  acting  on  behalf  of  the 
Commission: 

A.  Makes  any  warranty  or  representation,  express  or  implied,  with  respect  to 
the  accuracy,  completeness,  or  usefulness  of  the  information  contained  in 
this  report,  or  that  the  use  of  any  information,  apparatus,  method,  or 
process  disclosed  in  this  report  may  not  infringe  privately  owned  rights;  or 

8.  Assumes  any  liabilities  wiirh  respect  to  the  use  of,  or  for  damages  resulting 
from  the  use  of  any  information,  apparatus,  method,  or  process  disclosed 
in  this  report. 

As  used  in  the  above,  "person  acting  on  behalf  of  the  Commission"  includes 
any  employee  or  contractor  of  the  Commission,  or  employee  of  such  contractor, 
to  the  extent  that  such  employee  or  contractor  of  the  Commission,  or  employee 
of  such  contractor  prepares,  disseminates,  or  provides  access  to,  any  information 
pursuant  to  his  employment  or  contract  with  the  Commission,  or  his  employment 
with  such  contractor. 


UNCLASSIFIED 


Magneto-Pluld  Dynamics  Division 
Courant  Institute  of  Mathematical  Sciences 
New  York  University 


TID-4500 

35rd  edition  Physics 

MP-^1  NYO-1480-7 

Asymptotic  Equivalence 
of  the  Navler-Stokes  and 
Nonlinear  Boltzmann  Equations 

harod:  GRAD 

September  30,  1964 

AEC  Research  and  Development  Report 


Preprint  of  a  paper  to  be  published  In  the 

Proceedings  of  the  American  Mathematical  Society 

Symposium  on  Application  of  Partial  Differential 

Equations  in  Mathematical  Physics, 

New  York  City,  April  19^4. 


Contract  No.  AT(30-1 ) -l480 


-  1  - 
UNCLASSIFIED 


Abstract 

Previous  theory  of  the  linear  Boltzmann  equation 
Is  extended  to  somewhat  nonlinear  problems.   For  suf- 
ficiently small  Initial  data,  solutions  are  shown  to 
exist  for  macroscoplcally  long  times.   If  the  mean 
free  path  Is  small  the  nonlinear  solution  Is  particu- 
larly well  approximated  by  the  linear  solution,  and  It 
Is  also  approximated  by  the  macroscopic  Navler-Stokes 
solution  of  a.n  equivalent  problem. 
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1 .   Introduction 

One  of  the  most  interesting  features  of  the   use  of 
Boltzmann's  equation  to  describe  the  behavior  of  a  gas  Is  the 
possibility  of  a  radically  different  macroscopic  description 
of  the  gas  via  the  Navler-Stokes  or  Inviscld  Euler  Equations. 
The  fact  that  there  is  an  overlapping  range  where  both 
theories  are  believed  to  be  valid  on  physical  grounds  suggests 
the  strictly  mathematical  question  of  a  connection  between  the 
theories  of  the  Boltzmann  and  Navler-Stokes  equations.   The 
relationship  is  extremely  singular,  first  because  the  variables 
used  to  describe  the  state  of  the  gas  are  so  dissimilar,  and 
second  because  the  time  scales  on  which  the  gas  evolves  are 
apparently  unrelated. 

To  exhibit  this  connection  turns  out  to  require  much 
more  accurate  estimates  of  solutions  than  those  given  by  any 
naive  proofs  of  existence.   There  are  strong  nonlinear  existence 
theorems  in  the  large  for  the  spatially  homogeneous  problem. 
There  is  one  by  Carleman  [1]  for  hard  spheres,  one  by  Wild  [2] 
modified  by  Morgenstern  [3]  for  a  cutoff  Maxwelllan  potential; 
and  there  is  an  existence  theorem,  for  an  artificial  equation 
which  reduces  to  a  legltlm.ate  Boltsm-ann  equation  (with  a  general 
cutoff  potential)  in  the  space -independent  case  by  Povzner  [4]. 
But  this  space -independent  problem,  although  mathematically 
substantial,  is  physically  trivial  since  there  is  no  change  at 
all  in  any  of  the  macroscopic  variables  and  there  is,  of  course. 
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no  fluid  dynamics  without  spatial  variation.   Previous  existence 
theorems  In  problems  with  space  dependence  are  relatively  poor. 
Existence  Is  shown  for  a  time  comparable  to  a  collision  time; 
this  Is  quite  Irrelevant  on  any  macroscopic  time  scale  In  which 
fluid  behavior  Is  to  be  observed.   These  limitations  In  the 
nonlinear  spatially-dependent  theory  are  not  surprising  when 
one  notes  that  this  theory  should  be  expected  to  be  at  least 
as  difficult  as  the  corresponding  nonlinear  theory  of  the 
Navler-Stokes  equations  (which  should  arise  as  a  specal  limiting 
case ) . 

There  are  statements  in  the  literature  that  one  should  not 
even  expect  solutions  in  the  large  for  the  nonlinear  Boltzmann 
equation,  and  more  tractable  abstract  replacements  have  there- 
fore been  suggested  for  the  Boltzmann  equation  [7]j[^]-   Of 
course,  there  is  no  a  priori  reason  to  expect  the  existence  of 
solutions  in  the  large  for  any  nonlinear  equation.   This  doubt 
holds  just  as  well  (in  the  absense  of  any  information)  for  the 
equations  of  fluid  dynamics.   But  it  is  exactly  the  connection 
between  the  Boltzmann  equation  and  fluid  dynamics,  especially 
the  Navler-Stokes  equations,  that  leads  one  to  expect  the 
likllhood  of  existence  of  nonlinear  solutions  of  the  Boltzmann 
equation  in  the  large  (at  least  for  sufficiently  small  mean 


1.   See  Ref.  [5],  Chap.  III.   The  posthumous  proof  of  Carleman 
(Ref.  [6],  pp.  87-92)  is  vitally  Incomplete. 
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free  path).   As  a  matter  of  fact,  the  presence  of  strong 
results  in  the  spatially  homogeneous  case  where  fluid  dynamics 
Is  absent  suggests  that  nonlinear  estimates  similar  to  those 
In  fluid  dynamics  are  exactly  what  Is  missing  from  the 
Boltzmann  theory. 

According  to  some  recent  usage,  the  term  "Boltzmann 
equation"  Is  taken  to  refer  to  any  equation  which  Involves  a 
distribution  function  In  velocity  as  well  as  physical  space. 
We  shall  be  more  restrictive  and  consider  the  original  equation 
of  Boltzmann  and  Maxwell  In  Its  application  to  a  classical  gas. 
Even  with  this  restriction,  we  deal  with  a  general  type  rather 
than  with  a  specific  equation.   A  common  property  of  all  these 
equations  Is  that  of  conservation  of  mass,  momentum,  and 
energy.   This  Is  violated,  for  example,  by  equations  governing 
a  neutron  gas.   The  type  of  result  which  Is  usually  adequate 
In  the  latter  case,  viz.,  exponential  decay  or  exponential 
growth,  requires  relatively  trivial  estlm.ates.   In  the  linear 
version,  the  domilnant  feature  Is  whether  the  collision  opera- 
tor Is  positive  or  not.   For  a  classical  gas  the  operator  Is 
non-negative;  but  the  origin  belongs  to  the  spectrum  as  a 
consequence  of  conservation.   In  the  spatially  homogeneous 
case  the  zero  eigenvalues  can  be  discarded  and  we  do  obtain  the 
uniform,  exponential  decay  of  all  solutions.   In  the  general 
case  we  find  (but  this  is  a  much  more  subtle  result)  an  approx- 
imately exponential  decay  to  a  situation  in  which  fluid  behavior 
takes  over.   The  state  of  the  gas  then  evolves  at  a  much  slower 
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rate  and  ultimately  does  decay  to  equilibrium  despite  the  zero 
eigenvalues. 

The  variety  of  true  Boltzmann  equations  Is  associated 
with  the  choice  of  Intermolecular  force  law.   There  Is  at  the 
present  time  no  existence  theory  whatsoever  for  any  equation 
with  a  potential  that  extends  to  Infinity  (infinite  total 
cross-section).   This  Is  even  true  for  a  Maxwelllan  potential 
for  which  the  collision  operator  Is  known  to  have  a  complete 
pure  point  spectrum  with  explicit  polynomial  elgenfunctlons 
[8].   The  one  exception  Is  the  spatially  homogeneous  problem 
In  which  the  solution  Is  explicit  as  a  sum  of  exponentially 
decaying  elgenfunctlons.   Nothing  Is  known  about  the  spectrum 
for  any  other  non-cutoff  potential,  although  the  spectrum  Is 
known   for  a  large  class  of  cutoff  potentials  (and  Is,  in 
general,  not  discrete )[ 9] •   Even  with  cutoff  potentials,  the 

mathematical  tractablllty  seem.s  to  be  extremely  sensitive  to 

2 
the  type  of  cutoff.    But  there  does  exist  a  simple  way  of 

truncating  the  collision  cross-section  which  has  recently  been 
found  to  be  amenable  to  theory  [9]-   The  question  whether  the 
Boltzmann  equation  for  molecules  with  infinite  total  cross- 
section  is  mathematically  sound  (even  in  the  linear  case)  is 
completely  open. 

Discouragement  with  the  Boltzmann  equation  has  led  to  the 
study  of  various  types  of  modified  equations.   One  type  of 


2.   Personal  communication  from  L.  Pinkelstei 
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modification,  guided  by  mathematical  considerations,  replaces 
the  collision  operator,  which  is  classically  an  operator  in 
velocity  space  alone,  by  an  operator  which  is  nonlocal  in 
physical  space.   This  allows  proof  of  nonlinear  existence 
theorems  in  the  large  [■4-],  [7]-   But  these  models  seem  to  make  any 
property  other  than  the  bare  fact  of  existence  even  more 
elusive  than  originally.   For  example,  the  identification  with 
physical  parameters  is  hidden;  the  collision  operator  would 
require  re -definition  near  a  boundary;  the  relation  to  fluid 
dynamics  is  not  simple  and  is,  in  any  case,  not  that  of  an 
ideal  gas.   And,  of  course,  the  basic  arbitrariness  in  the 
kernel  describing  the  space -averaging  is  awkward.   If  the  solu- 
tions could  be  shown  to  be  insensitive  to  the  precise  degree 
of  averaging  in  physical  space,  then  a  theory  of  the  actual 
Boltzmann  equation  would  probably  ensue.   This  type  of  uniform 
estimate  would  seem  to  be  comparable  in  difficulty  to  the 
uniform  estimates  which  are  required  to  establish  the  singular 
connection  with  fluid  dynamics.   At  present  there  are  no  such 
uniform  estimates,  and  the  future  status  of  these  models  is 
uncertain. 

An  entirely  different  type  of  modification  of  the  Boltzmann 
equation,  referred  to  as  a  relaxation  model,  is  guided  by  physi- 
cal principles  and  an  entirely  different  type  of  mathematical 
expediency  (e.g.  see  [10]).   These  equations  exhibit  no  greater 
simplicity  with  regard  to  general  existence  theory,  but  do  offer 
practical  advantages  of  explicit  solutions  in  simple  problems 


and  are  more  attractive  for  numerical  work. 

Existence  of  solutions  to  the  linear  Boltzmann  equation 
for  hard  spheres  was  shown  by  Carleman  (Ref.  [6],  Part  2, 
Chap.  I).   Although  existence  was  proved  In  the  large,  the 
growth  estimate  e   ' '^    (t  Is  the  mean  collision  time)  Is  cata- 
strophlcally  poor  on  a  macroscopic  time  scale  (typically 
T  ^v_  10~   sec.).   This  linear  existence  theory  has  been 
generalized  recently  to  Include  a  wide  class  of  Intermolecular 
potentials,  and  It  has  been  Improved  to  show  boundedness  for 
all  time  [9]  and  even  an  ultimate  decay  to  equilibrium  [11]. 
Even  more  Important,  the  bounds  have  been  shown  to  be  Indepen- 
dent of  the  mean  free  path  (which  enters  as  a  singular 
parameter)  and  therefore  carry  over  to  the  macroscopic  regime. 
These  strong  estimates  have  been  used  to  show  the  relation 
between  Boltzmann  solutions  and  fluid  solutions,  making 
rigorous  the  Hllbert  theory  [12];  and  In  this  paper  we  use 
the  strong  linear  estimates  to  derive  certain  nonlinear  exis- 
tence results  and  to  evaluate  the  legitimacy  of  the  linear 
theory. 

The  connection  between  the  Boltzm.ann  equation  and  fluid 
dynamics  was  shown  formally  by  Hllbert  [13].-^   Introducing  a 
singular  parameter  e  Into  the  Boltzmann  equation  (e  is  the  mean 
free  path  or  mean  collision  time,  cf.  Sect.  2),  Hllbert  showed 
that  an  expansion  of  the  solution  in  powers  of  e  has  the  formal 


3-  For  a  critical  account  of  the  Hllbert  and  Chapman -En skog 
theories  see  Ref. [5],  Chap. IV;  for  recent  developments  see 
Ref  .[12]  . 
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property  that  each  term  in  the  expansion  Is  determined  by  an 
algorithm  Involving  the  fluid  state  alone.   A  plausible  argu- 
ment can  be  given  to  show  that  Hllbert's  expansion  Is  asymptotic 
In  e  except  for  a  short  Initial  layer,  a  narrow  boundary  layer, 
and  any  Internal  shock  layers  (Ref.  [5],  ^  26).   The  linear 
estimates  mentioned  above  (uniform  In  s)  provide  a  proof  that 
the  Hllbert  expansion,  suitably  modified.  Is  asymptotic  to 
arbitrary  solutions  of  the  linear  Boltzmann  equation  with  smooth 
Initial  data  [12].   Specifically,  after  an  Initial  transient  of 
duration  0(e),  the  distribution  function  becomes  approximately 
locally  Maxwelllan.   The  fluid  state  is  governed  by  the  Invlscld 
Euler  equations  of  fluid  dynamics  for  a  time  0(l)  and  by  the 
Navler-Stokes  equations  for  a  time  0(l/e)  (ths  viscosity  coef- 
ficient Is  0(e)). 

In  this  paper  we  use  the  linear  estimates  to  prove  existence 
of  solutions  to  the  nonlinear  Boltzmann  equation  for  a  length  of 
time  which  depends  on  the  amplitude  of  the  Initial  data.   For  suf- 
ficiently small  Initial  amplitude,  the  solution  exists  for  Inter- 
esting macroscopic  times,  but  It  Is  In  general  approximated  by  the 
linear  equation  for  only  a  short  (microscopic)  part  of  this  time. 
If  the  mean  free  path  Is  small  In  addition,  then  the  linear  approx- 
imation Is  valid  for  much  longer  macroscopic  times.   Furthermore, 
both  solutions  are  approximated  by  fluid  dynamics  (Invlscld 
Euler  or  Navler-Stokes  equations  depending  on  the  circumstances). 
The  common  fluid  behavior  for  small  e  Is  the  reason  for  the 
Increased  range  of  validity  of  the  linear  approximation.   The 
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natural  time  of  evolution  of  fluid  properties  is  macroscopic 
and  bears  no  relation  to  the  mean  collision  time,  which  is 
the  basic  time  scale  for  Boltzmann  solutions. 

The  principal  mathematical  difficulty  in  the  linear  theory 
arises  from  the  fact  that  the  linear  collision  operator, 
although  non-negative,  is  not  positive  definite.   But  it  is 
exactly  this  property  which  provides  the  fluid  behavior  in 
the  limit  of  small  e.   There  is  a  rapid  approach  in  a  time  of 
order  e  toward  a  state  which  is  essentially  macroscopic.   This 
is  followed  by  fluid  behavior,  largely  inviscid  on  a  finite 
time  scale,  but  viscously  decaying  on  a  long  time  scale  of 
order  l/e.   Without  the  zero  eigenvalues,  all  solutions  would 
decay  to  zero  in  a  time  comparable  to  e  and  both  the  theory  and 
the  physical  behavior  would  be  relatively  trivial. 

The  principal  difficulty  in  the  nonlinear  problem  is  that 
a  priori  estimates  in  the  linear  problem  are  natural  in  a  norm 
which  is  inappropriate  for  estimation  of  the  nonlinear  terms. 
To  find  estimates  in  compatible  norms  requires  the  use  of 
explicit  smoothing  properties  of  a  part  of  the  linear  collision 
operator  (after  removing  a  singular  and  unbounded  part)  and 
depends  heavily  on  the  assumed  rectangular  domain  and  simple 
boundary  conditions.   The  rectangular  domain  with  perfectly 
reflecting  walls  has  the  special  property  of  eliminating  the 
boundary  layer  which  was  mentioned  earlier  and  which  requires 
special  estimates.   The  boundary  layer  is  of  course  absent  in 
an  unbounded  domain.   But  this  case  requires  additional  estimates 
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of  the  decay  at  Infinity  In  physical  space;  this  problem  Is 
treated  elsewhere. 


5 
2 .   Properties  of  the  Boltzmann  Equation 

We  consider  a  dlmenslonless  distribution  function  F(^,x,t) 
of  the  dlmenslonless  arguments  velocity  (^),  space(x),  and  time 
(t).   The  variables  have  been  made  dlmenslonless  with  respect 
to  the  parameters  of  a  reference  state  which  Is  specified  by  a 
mass  density  p   and  mean  thermal  speed  C   =  (RT  )^.   The  dis- 
tribution function  has  been  made  dlmenslonless  by  a  factor 
p  /C  ,  the  arguments  ^,   x,    t  by  respective  factors  C  ,  a  mean 

free  path  A   (related  to  p  ,  C   and  molecular  constants),  and 
^      o  ^  "^0    0  ^  ' 

a  mean  free  time  ^^VC  .   A  skeleton  form  of  the  Boltzmann  equa- 
tion is 

The  collision  term  Q,  is  a  quadratic  operator  on  P  as  a  function 
of  I  alone;  the  streaming  term.  DF  operates  on  F  as  a  function 


4.  To  appear  in  Communications  on  Pure  and  Applied  Mathem.atlcs . 

5.  For  a  more  detailed  description  and  proofs  of  the  properties 
quoted  in  this  section  see  Ref.  [9]- 
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of  (x,t)  alone  with  ^  as  a  parameter.   Initial  and  boundary 
values  can  be  expected  to  be  the  same  as  for  the  differential 
operator  alone. 

The  dimenslonless  equilibrium  solution  (Maxwellian)  which 
corresponds  to  the  reference  state  (p  ,  C  )  is  F  =  co(|)  where 


l^^exp(-i|2^   .  (2.2) 

(27r)- 


CD  =  yp^    tixp^,-  2 


We  have 

Q(co,a3)  =:  0   .  (2.5) 

We  are  interested  in  solutions  of  (2.1)  which  are  not  greatly 

perturbed  from  F  =  co;  to  this  end  we  introduce 

\_ 
F  =  CO  +  0)2  f  (2.4) 

\_ 
(the  factor  co^  is  convenient  for  normalization  of  f).   The 

equation  staisfied  by  f  is 

DF  +  LF  =  vn(f,f)  (2.5) 


where 

LF  =  -  2co"2  Q(a3,aj2  f)  (2.6) 


i       1 

"2 


and  the  function  v(^)  will  be  defined  in  a  moment.   The  linear 
operator  L  is  singular  and  unbounded.   However,  for  a  large 
class  of  appropriately  restricted  intermolecular  potentials,  it 
can  be  shown  to  split  into  two  terms, 
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Lf  =  vf  -  Kf  .  (2.8) 

The  singular  and  unbounded  part  of  L  Is  given  by  the  first  terra; 
K  Is  a  bounded  (even  compact)  Integral  operator.   The  factor 
v(^),  the  collision  frequency  of  a  molecule  of  speed  I ,  Is  a 
certain  unbounded  positive  function  of  ^.   The  factor  v  in 
(2.7)  was  removed  because  the  remaining  factor  |  'will  be 
found  to  be  bounded  In  an  appropriate  sense. 

We  now  quote  a  number  of  properties  of  the  operator  L 
which  are  proved  elsewhere  [9];    the  required  properties  of  P  are 
derived  here  in  the  Appendix.   The  derivation  of  these  properties 
is  dependent  on  specific  assumptions  with  regard  to  the  inter- 
molecular  force  law;  mainly  that  it  be  repulsive,  of  finite  total 
cross-section  (suitably  defined),  and  sufficiently  "hard"  (i.e. 
the  force  is  sufficiently  strong  for  small  distances  of  separa- 
tion) .   Classical  hard  spheres  and  suitably  cut  off  power  law 
forces  with  exponent  >  5  are  Included. 

The  operators  L  and  K  are  symmetric  and  L  is  non-negative; 
these  statements  are  made  relative  to  the  Lp  norm  in  the 
infinite  three-dimensional  space  ^, 

llf  11^  =  /f^d^ 

(2.9) 
(f,g)  =  /  f  gd^  • 

A  critical  point  in  all  the  analysis  is  that  L  is  not  positive. 
The  origin  is  a  multiple  eigenvalue 

L^^  =0    (a  =  0,1,2,3,^)  (2.11) 
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with  the  explicit  normalized  elgenfunctlons 


i/    =00 


o 
if^   =   ^^03*    (1  =  1,2,3)  (2.11) 

\  =    {(^   -3)(co/6)*   . 

The  existence  of  these  elgenfunctlons  Is  closely  related  to  the 
fact  that  mass,  momentum,  and  energy  are  conserved;  and  It  Is 
because  L  Is  not  positive  that  solutions  of  the  Boltzmann  equa- 
tion do  not  decay  simply  to  zero  but  exhibit  fluidlike  behavior 
under  appropriate  circumstances. 

The  function  v(£,)  is  usually  monotone,  it  is  bounded  away 
from  zero  (for  the  class  of  hard  potentials  under  consideration), 
and  it  is  usually  unbounded  for  large  |  as  a  power  I ',  0  <  y  <  1. 
We  shall  require  only  the  bounds 

1  <  v(^)  <  v^(l+|2)*  (2.12) 

where  v   Is  a  constant.   The  lower  bound  has  been  taken  to  be  1 
by  a  suitable  choice  of  the  normalizing  dimensional  mean  collision 
time . 

Since  K  is  compact,  the  essential  spectrum  of  L  consists  of 
exactly  the  values  taken  by  v(^).   Therefore  the  origin  is  an 
Isolated  point  of  the  spectrum,  and  we  have 

(f,  Lf)  >  ^L(f,f) 

(2.13) 
(f,  V^J  -   0 
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for  some  positive  constant  [i. 


We  denote  by  k  the  bound  of  K, 


IlKf  II  <  K^  ||f  II  .  (2.14) 

Other  more  refined  properties  of  K  are 

max  |Kf  1  <  K-^  IJf  II  (2.15) 

max(l+^^)^''+^)/^|Kf|  <  K^^2"^ax(l+^2)''/2|f|       (2.l6) 
(r  =  0,  1,  ...)• 


and 


In  other  words,  Kf  Is  bounded  If  f  Is  square  Integrable,  and  Kf 
decays  faster  than  f  by  a  factor  £,  for  large  ^. 

There  are  several  alternate  forms  of  the  Boltzmann  equation 
that  we  shall  use.   If,  instead  of  (2.4),  we  introduce 

P  =  CO  +  ecD^f  (2.17) 

into  the  Boltzmann  equation  (2.1),  we  obtain 

Df  +  Lf  =  evr'(f,f)  .  (2.18) 

If  e   is  small,  we  are  led  to  consider  the  linear  equation 

Df  +  Lf  =  0  (2.19) 

with  the  expectation  that  solutions  of  this  equation  will  approx- 
imate those  of  (2.18)  for  small  e. 

A  more  interesting  change  of  scale  is 

P  -  cD/e  +  oo^f  (2.20) 
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which  leads  to  the  equation 

Df  +  I  Lf  =  vr(f.f)  •  (2-21) 

Implicit  in  (2.20)  is  the  supposition  that  the  perturbation  is 
not  only  small  relative  to  the  reference  state  co/e,  but  that 
the  perturbation  has  a  basic  length  scale  which  is  long  compared 
to  the  mean  free  path  of  the  reference  state  [which  is  0(e) 
since  the  density  is  0(l/e)].   It  is  apparent  that  e  =  0  is  a 
singularity.   Moreover,  it  is  not  Immediately  evident  that  solu- 
tions of  (2.21)  will  be  approximated  by  those  of  the  associated 
linear  equation 

Df  +  -  Lf  -  0  .  (2.22) 

£ 

Despite  appearances,  it  will  turn  out  that  solutions  of  the 
linear  equation  (2.22)  do  approximate  those  of  the  nonlinear 
equation  (2.21).   The  reason  is  their  common  fluidlike  behavior 
in  the  singular  neighborhood  of  small  e. 

A  final  change  of  scale  which  we  shall  wish  to  consider  is 

i_ 

which  leads  to  the  equation 

Df  +  -  Lf  =  evp(f,f)  .  (2.24) 

This  equation  is  "linearized"  to  a  higher  order  (e  )  than  the 
mean  free  path  is  shortened  (e).   It  is  associated  with  the  same 
linear  equation  as  (2.21),  viz.  (2.22),  but  the  relation  is,  in 
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this  case,  more  evident. 

Implicit  in  the  various  scalings  in  e  is  the  assumption 
that  the  initial  data  are  independent  of  e. 


3.   Boundary  Conditions 

We  consider  a  rectangular  domain  and  two  types  of  boundary 
conditions,  periodicity  and  specular  reflection.   The  simpler 
case  is  periodicity.   We  take  a  fundamental  domain  with  sides 
(a-,,  ap,  a^)  and  assume  that  the  initial  values  are  periodic 
in  X.  with  periods  a. .   It  is  easily  verified  that  this  peri- 
odicity is  retained  by  all  functions  that  arise  in  the  analysis, 
in  particular  by  the  solution.   The  periodicity  property  is  used 
essentially  at  only  one  point  in  the  analysis,  to  show  that  a 
boundary  integral  vanishes  in  m.aking  an  a  priori  estimate. 

The  specular  reflection  boundary  condition,  also  applied 
in  a  rectangular  domain,  states  that  f  is  an  even  function  of 
the  normal  com.ponent  of  velocity  at  a  wall 

f(l)  =  f(^  -  2n(^n))  .  (3.1) 

By  reflection  of  the  fundamental  domain  D  with  respect  to  each 
of  three  coordinate  planes  we  obtain  a  domain  D*  consisting  of 
eight  replicas  of  D.   In  D*  the  function  f  satisfies  a  perio- 


18 


dlclty  boundary  condition.   It  Is  easy  to  verify  that  the  sym- 
metry of  the  Initial  function  f  with  respect  to  the  three  planes 
which  bisect  D*  is  maintained  by  all  functions  which  arise  in 
the  analysis.   Thus  solution  of  the  Boltzmann  equation  with  a 
periodic  boundary  condition  in  D*  yields  the  desired  solution 
with  specular  boundary  condition  in  the  subdomain  D. 

There  are  a  few  precautions  that  should  be  observed  in  passing 
between  D  and  D*.   If  the  Initial  function  f  is  required  to  be 
continuous  in  D*,  this  implies  that  the  boundary  condition  (3-1) 
must  be  satisfied  by  the  initial  function  in  D.   A  continuous 
initial  derivative  Sf/3x-,  in  D*  implies  the  condition 


dxl 


f(^l)  =  -^f(-^l)  (5-2) 


at  the  wall  perpendicular  to  x-,  ,  etc.   In  other  words,  smoothness 
in  D  Involves  satisfaction  of  certain  relations  beyond  conven- 
tional smoothness  at  the  boundary.   These  relations  can  be 
Interpreted  as  smoothness  conditions  as  seen  by  a  given  particle 
which  is  specularly  reflected  (and  changes  its  velocity  dlscon- 
tlnuously  while  doing  so).      : 

The  problems  of  initial  data  which  do  not  satisfy  (3.1)  or 
(3-2),  of  more  general  discontinuous  initial  data,  and  of  general 
nonrectangular  dom.ains  are  closely  related  but  will  not  be  con- 
sidered here. 
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k .      Norms 

In  order  to  properly  combine  estimates  for  the  linear  and 
nonlinear  equations  It  Is  necessary  to  Introcude  an  assortment 


of  norms.   We  denote  by  N  the  Lp  norm  In  (l,x)  (as  distinguished 


from 


,  previously  used  for  the  norm  In  2.) 


N  [f]  -  ■^  /  f^d^  dxp 


(4.1) 


The  simple  maximum  norm  we  denote  by 


N  [ f ]  =  max  I f 
o       t  ' 

e,  ,x 


(4.2) 


and  more  generally 


N  [f]  =  max  (l+^^)'^/^|f 
l.x 


(^•3) 


We  shall  also  use  the  combined  norms 


*f^[f]  =   max  (1+1^)^/^1  /f^dxl^ 


(4.4) 


and  the  norms  of  x-derlvatlves 


N^[f]  =  max  N  [\/^f] 
n=o 

N^[f]  =  mix  N^[V"f] 
^      n=o 

w![f  ]  =  mix  ?^-^[V->  ] 
n=o 


(4.6) 


where 


V"f|2  .  2 


1-,  .  .  .  1 
1    n 


S^f 


dx.     77".    dx. 

^1 


(4.8) 
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We 


shall  use  the  symbol  "^  to  represent  any  one  of  the  norms 


.s   ,,s 


N,  N  ,  ¥  .   Each  of  the  norms  Is  evaluated  at  a  given  time,  t 
For  the  maximum  with  respect  to  t  we  use 


1q=   max  ^(t)  (4.9) 

where  max  is  understood  to  be  taken  for  0  <  t  <  t  ,  t   finite 
or  infinite  depending  on  the  context. 
Some  elementary  estimates  are 

¥^  <  w  ^, ,     ^l  <  nJ"^^  (4-11) 

r  —  r+1        r  —  r  ^  ' 

N  <  TT^g  (4-12) 

2  -1 
the  last  being  a  consequence  of  the  fact  that  (l+l  )    is  square 

2-2       2 
integrable   [  /  (1+^  )   d|  =  tt  ]  .    There  exists  a  constant  s 

depending  on  the  domain  such  that  (Sobolev) 

N^  <  s  *f^+2.  (4-13) 

r    o   r  \        ^  / 

Prom  Sect.  2  we  transfer  the  Lp  inequalities  (2.13)  and 
(2.l4)  unaltered  to  the  space  (^ ,x)  and  generalize  (2.15)  and 
(2.l6)  to  read 

¥^[Kf]  <  K^N^[f]  (4 -14) 

¥^^^[Kf]  <  K^+2*^^[f]  .  (4.15) 


-  21  - 


Similar  estimates  for  the  quadratic  collision  term  \~' , 
derived  In  the  Appendix,  are 


and 


N^[  r(f  ,g)]  <  YQN^[f]N^[g]  (4.16) 


*f;^[r(f,g)]  <  YoW3.[f]N^[g]  .  (4.17) 


There  Is  no  estimate  of  [  In  terms  of  only  Lp  estimates  of 
f  and  g. 

The  above  estimates  of  K  and  |  '  form  the  basis  for  the 
entire  theory.  We  continue  with  some  consequences  of  these 
estimates . 

From  K'Vfl  <  (1+^  )^IVf|  we  obtain 

*^J[^Vf]  <  Nj+^[f]  (4.18) 

and  from  (2.12)  we  have 

*fj[vf]  <  v-^*fj^^[f]  (4.19) 

from  which  follows 

$f^[Lf]  <  (v^  +  K^+i)<.+i[f]  (4.20) 

using  (2.8),  (4.15),  and  (4.11).   If  f  Is  a  sufficiently  regular 
solution  of  the  linear  Boltzmann  equation  (2.19),  Sf/5t  can  be 
estimated  from  (4.l8)  and  (4.20)  as 

W^[af/St]  <  (1  +  v^  +  ^r+l)*^Jti^^^  •  ^^-2^) 

A  vital  property  of  all  norms  is 
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riCf]  <  7i[f]  (^-22) 

where 

t{i,x,t)   =  f(^,x-se,t)  (4.23) 

and  s  is  a  constant  (note  that  this  would  not  be  true  if  ^  were 
replaced  by  the  stronger  norm  with  max  preceding  /•■•  dx). 

The  basic  formula  to  which  all  convergence  estimates  will 
be  referred  is  the 
Lemma :   If 

f(^,x,t)  =  ^   e-v(t-s)/e  g(^^x-|(t-s),s)ds      (4.24) 


then 


e 
6 


^r^^'^   ^  ^r^^s]  •  (4.25) 


If  V  is  a  constant,  this  estimate  for  any  norm  is  classical 
since  f  is  a  mean  value  with  parameter  s  and  total  weight  less 
than  one  of  the  functions  g, 

t 


(v/e) 


r 


e 


-v(t-s)/e  ^  -^  _  g-vt/e  ^  ^ 


o 


With  V  a  function  of  ^,    the  estimate  (4.25)  is  not  true  for  the 
L2  norm  of  N  and  a  maximum  norm  with  respect  to  i    is  needed. 
For  the  proof,  we  use  Schwartz's  inequality  to  obtain 

/  g(^,x-^(t-s^),s^)g(t,x-^(t-S2)  ,S2)dx 

<  (i+^^)-^(i^[g])2 

which  implies 

/  f2dx<  {l+er^{i^[g]f{l-e-'''/^f    , 
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and  thus  (4.25) • 


5-   Outline  of  the  Proof 

The  linear  equation  (2.22)  can  be  written 

11+  i-Vf   +  ^  vf  =  1  Kf  (5.1) 

and  transformed  into  the  Integral  equation 

f(^,x,t)  =  f^(^,x-et)exp(-vt/e) 


1 

+  — 

e 

o 


exp[ -v(t-s)/e]k(^,x-4(t-s),s)ds       (5-2) 


where 


and 


k  =  Kf  (5-3) 


f(|,x,0)  =  f^(^,x)  (5.4) 

Is  a  given  initial  function.   We  solve  (5-1)  by  iteration  on  the 
right  side,  or  the  equivalent.  Iterate  (5-2)  and  (5-3)  succes- 
sively.  The  crude  estimate  obtained  by  setting  the  exponentials 
in  (5.2)  equal  to  one  yields  convergent  iterations  with  a  growth 
exp(Kt/£)  in  any  norm  /?  where  k  is  the  bound  of  K  in  the 
particular  norm  ?] .      No  great  improvement  in  this  growth  esti- 
mate can  be  made  by  examining  the  iterates  more  carefully 
(unless  V  is  a  constant,  see  below)  because  of  the  nontrivial 
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nullspace  of  the  operator  L;  Kf  is  "just  as  large"  as  vf. 

But  for  any  smooth  solution  f ,  we  can  obtain  the  a  priori 
estimate 

N[f]  <  N[f^]  .  (5.5) 

Multiplying  through  the  equation  (5-1)  by  f  and  integrating 
with  respect  to  x  and  i,   we  can  drop  the  streaming  term 
V"(  2^f^)  using  periodicity  in  x  and  estimate  the  collisional 
'  term  by  (f,  Lf)  >  0,  to  obtain  (5-5) •   A  simple  argument  shows 
that  (5.5)  obtains  without  unnecessary  smoothness.   In  only 
the  case  v  =  constant  (cutoff  Maxwellian  molecules)  can  the 
uniform  bound  on  N[f]  be  established  directly  for  each  itera- 
tion (the  Lp  norm  of  K  is  just  v).   If  v  is  a  function  of  i, 
we  are  blocked  by  the  inapplicability  of  the  Lemma,  (4.25),  in 
an  Lp  norm. 

The  mollifying  properties  of  K,  (2.15)  and  (2.l6),  can  now 
be  used  to  show  boundedness  of  the  solution  as  a  function  of  ^, 
and  a  Sobolev  lemma  to  show  boundedness  as  a  function  of  x  for 
sufficiently  smooth  data  (all  x-derivatlves  of  f  satisfy  the 
same  equation  and  boundary  conditions  and  therefore  the  same 
a  priori  estimate).   Note  that  all  these  unlf orm-in-tlme  esti- 
mates are  also  independent  of  the  value  of  e. 

Similar  estimates,  also  Independent  of  e,  can  be  found  for 
the  Inhomogeneous  linear  equation.   Using  them,  we  iterate  on 
the  right  side  of  the  nonlinear  equation  (2.21).   The  unbounded 


6.  This  is  essentially  the  linearized  version  of  Boltzraann's 
H-theorem. 
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factor  V  on  the  right  Is  taken  care  of  by  the  strong  Inhomo- 
geneous  linear  estimate  which  gives  a  solution  decaying  faster 
In  ?,  than  the  Inhomogeneous  term.   A  more  serious  problem  is  to 
obtain  maxlm.um  estimates  which,  for  the  linear  equation,  are 
obtained  from  the  a  priori  Lp  estimates  only  after  losing  derlva- 
tlves.   This  Is  taken  care  of  by  using  a  norm  ^^.   Writing 
f.  =  bf/bx    ,      f   =  3  f/bx.bx . ,    etc.,  we  have 

g|^r(f,f)  =  2r(f,fi) 

^^  r(f,r)  =  2r{f,fij)  +  2r{fi,fj)  (5.6) 

I     1   J   k  ^  ^ 

i  +  2p(f  .,f,  .  )  +  2p(f,  ,f .  .)  . 

'  ^  J '  ki^     '  ^  k   ij  ^ 

An  estimate  for  w[f],   using  Sobolev,  Implies  a  maximum  esti- 
mate N_[f].   Each  term  on  the  right  side  In  (5-6)  has  at  least 
one  argument  of  ]    'which  Is  f  Itself  or  a  first  derivative;  i.e. 
at  least  one  argument  is  bounded,  and  the  other  one  is  Lp .   Thus 
t4-  [  l~'(f  ,f )  ]  can  be  estimated  in  terms  of  ^^[f].   In  this  way 
existence  is  proved  for  (2.18) 

11+  e-Vf  +  Lf  =  £vr(f,f)  (5.7) 

for  a  time  of  order  l/e  and  for  (2.21) 

U+  l-Vf  +  ^  Lf  ==  vr(f.x)  (5.8) 

for  a  time  of  order  1.   The  two  results  are  equivalent  in  that 
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each  describes  existence  for  a  long  time  of  order  l/e  times 
the  mean  collision  time.   But  the  class  of  initial  data  in 
(5.8)  is  more  special  in  that  the  f^  is  not  only  small  but  it 
is  also  slowly  varying  in  space  compared  to  the  mean  free 
path.   We  may  obtain  (5-8)  directly  from  (5-7)  by  a  change  of 
scale,  viz.,    by  inserting  f(ex)  and  changing  the  time  scale 
to  et .   In  return  for  the  greater  specialization  of  initial 
values  in  (5.8),  we  find  that  the  solution  is  more  accurately 
linear.   In  (5.7)  the  solution  is  approximated  by  the  corres- 
ponding linear  solution  (2.19)   to  within  an  error  0(e) 
microscopic  time  0(1);  this  is  evident  by  inspection  of  the  non- 
linear equation.   But  a  more  powerful  estimate  of  the  inhomo- 
geneous  linear  equation  can  be  made  in  (5-8)  and  yields  the 
unexpected  result  that  the  solution  of  (5-8)  is  approximated 
by  the  corresponding  linear  solution  (2.22)  for  the  entire 
(macroscopic)  time  of  existence.   The  proof  of  this  result  is 
closely  related  to  the  proof  that  the  solution  of  (5.8)  is 
approximated  by  fluid  dynamics.   Since  the  parameter  e  does  not 
appear  in  the  invlscid  fluid  equations, it  is  not  surprising  that  a 
small  fluid  perturbation  of  relative  size  e  will  remain  accurately 
linear  for  the  macroscopic  time  0(l).   For  the  "more  linear" 
equation  (2.2^),  existence  and  approximate  linearity  as  well  as 
fluid  behaviorare  proved  for  a  time  0(l/e).   The  comparisons 
with  fluid  dynamics  are  contained  in  an  earlier  paper  [12],  but 
we  shall  sketch  the  proofs  here  for  completeness. 

More  efficient  use  of  the  same  estim.ates  that  show  that 
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the  solutions  of  (5.8)  are  approximately  linear  would  give   an 
Improved  existence  proof  for  (5-8)  for  the  longer  time  O(log  1/ e) , 
but  we  shall  not  bother  with  this. 


6.   The  Linear  Equation ' 

Consider  first  the  Integral  form.ulatlon  (5.2),  (5-3)  with 
an  Initial  function  f  (^,x)  which  Is  periodic  and  bounded  with 


respect  to  one  of  the  norms  /(  .   Iterating  and  estimating  |f| 
by  setting  v  =  0, 

t 
^[f"+l]  <?/[fJ  +|J?/[k-(s)]ds 

(6.1) 


■.n 


we  obtain  a  contracting  sequence  of  Iterates  exp(-Kt/e)f'  In 
the  norm  77  where  k  Is  the  bound  on  the  operator  K  (k  =  k  for 
'}!  =.   N^,  K  =  K2_L  .  for  7?  =  N^  or  *f^)  .   We  may  choose  to  look 
at  (5.1)  with  V  =  0  as  a  dominating  equation. 

The  function  f(^,x,t)  defined  by  the  contracting  sequence 
Is  In  the  space  given  by  the  completion  of  /J   and  satisfies 
(5.2),  (5.3).   This  solution  is  unique  in  the  class  of  all 


7.  The  theory  of  this  section  is  taken  from  Ref.[9].   We  repeat 
the  proofs  (in  a  slightly  different  form)  to  make  the  abstract 
existence  theory  self  contained  (i.e.  starting  from  stated 
properties  of  the  collision  operators). 
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functions  f(|,x,t)  which  have  a  finite  norm  7?  for  some  inter- 
val 0  <  t  <  t  ;  and  at  any  finite  t,    there  is  continuous 
dependence  of  f  on  f   in  the  norm  /\  . 

The  solution  has  the  same  houndedness  and  differentiability 
properties  inherent  in  the  norm  TJ   as  the  initial  function. 
Let  us  first  assume  that  N[f  ]  is  finite.   From  (5-2)  we  deduce 
the  existence  of  the  directional  derivative  Df,  and  in  this 
generalized  sense  f  satisfies  the  original  equation  (5.I).   If 
^-.[f  ]  is  finite,  then  Sf/Bt,  t,-\/f,    and  vf  exist  individually 
in  (5.1)  (as  measurable  Lp  functions).   If  the  maximum  norm 

N  [f  ]  is  finite  and  f  as  well  as  the  indicated  derivatives 
r   o  o 

are  continuous,  then  the  solution  is  continuous  together  with 

its  derivatives  by  the  uniform  convergence  in  the  maximum  norm. 

Now  consider  a  fixed  initial  function  f  with  finite  norm 

o 

in) 
N.   We  can  approximate  f  by  a  sequence  of  functions  f^  '   with 

continuous  first  derivatives  and  finite  (but  not  uniformly 
bounded)  N^  [f^*^^]  such  that  N[f^-  f^^^]  — >  0.   By  the  continu- 
ous  dependence  on  initial  data  N[f^  ']   — >N[f]  at  any  time  t. 
But  each  of  the  smooth  solutions  f ^"''  satisfies  the  a  priori 
estimate  N[f^"'']  <  NLfJ^"-']-   This  follows  by  multiplication 
through  the  Boltzmann  equation  (which  is  satisfied  pointwise) 
by  f^^' ,    noting  that  ^-  V(f^^  )^  is  integrable  in  ^  and  con- 
tinuous.  Consequently  f  itself  satisfies  the  estimate 

N[f]  <  N[f^]  .  (6.2) 

Exactly  the  same  argument  holds  for  any  x-derivatives  of  f 
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which  are  square  integrable, 

N^[f]  <  N^[f^]  .  (6.3) 

For  the  linear  theory  alone  (cf.  Ref.  [9]),  It  Is  not  necessary 
to  be  so  parsimonious  of  derivatives  In  deriving  the  a  priori 
estimate;  but  It  Is  vital  for  the  nonlinear  application. 

To  extend  these  uniform  Lp  estimates  to  maximum  estimates, 
we  turn  to  the  Integral  equation  (5.2).  Referring  to  the  lemma 
of  Sect.  ^   and  recalling  that  v  >  1,  from  (5-2)  we  have 

ii-^.    [f]  <  Wj  [f^]  +  i^  [k]  .  (6.4) 

Using  (4.14),  (4.15)  In  (50),  we  find 

W^  [k]  <  K^N^[f] 

(6.5) 

^-  M  <   Vi*^5_i[f]  . 

and  Inserting  these  In  (6.4)  we  have 

^l     [f]  <  ^^  [f^]  +  K^F[f] 

(6.6) 

Solving  these  recurrence  Inequalities  for  successive  J  and 
Inserting  the  a  priori  estimate  (6.3), 

ill    ^^^  ^^1  ^^^o^  +  ^2*^0  f^o^  ^-   ^2^1^^t^o] 
<  (l+Kg)^^  [f^]  +  K2K^N^[fQ] 
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n|     [f]     <      (1    +     K^    +     K^K2)*f|     [f^]     +     K^K^K^N^Lf^] 
<     (1    +     K-^    +     K^Kg    +    TTK    K^^^)ii^     [f^] 

and   for  r  >   2 

where  a  depends  on  k,  .  .  .  k  ,  .,  only.   These  estimates  are 
r  1      r+i 

distinguished  from  the  JLg  estimates  {6.J>)    In  that  the  constants 
a  are  necessarily  larger  than  one. 
For  the  Inhomogeneous  equation 

Df  +  I  vf  =  I  Kf  +  vg  ,    ^o  "  °  ^^'^^ 

we  proceed  similarly.   Iterating 

,t 


f(l,x,t)  = 


o 


exp[ -v(t-s)/e]  ^  (^  ,x-|  (t-s)  ,s)ds 

(6.9) 


5  =  I  Kf  +  vg 

we  obtain  existence  of  f  over  any  time  Interval  In  which  the  norm 
of  vg  Is  finite.   As  an  a  priori  estimate,  using  the  fact  that 
L  Is  non-negative,  we  have 


r^o 


N^[f]  < 


J. 


N^[vg]dt  <  t^N^[vg]  <  TTV^t^^f^  [g]  .    (6.10) 


o 


We  obtain  maximum  estimates  from  the  Integral  expression  (6.9) 
by  use  of  the  Lemma  (4.25) 

*f^  [f]  <  K^N^[f]  +  ei^  [g] 

(6.11) 
^5  [f]  <  ViN^^_^[f]  ^e^'.    [g] 
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Iterating  these  Inequalities  and  taking  r  >  5,  we  find 

*fj  [f ]  <  (b^t^  +  £C^)J^J  [g]  ,    r  >  5  .  (6.12) 

Summarizing  the  essential  results  we  have 
Theorem  1:   There  is  a  unique  solution  to 

Df  +  -  Lf  =  vg 

for  ¥  [f  ]  and  #  [g]  finite  and  r  >  3  which  satisfies 

<  ^n   <  a^^^  [f^]  +  (b^t^  +  ec^)^l   [g]  .        (6  13) 

The  constants  a  ,  b  ;,  c   depend  only  on  v.,  ,  k,  .  .  .  k   ,  ,  and, 

rrr'^       o'     1^1     r+1 

in  particular,  are  Independent  of  e  and  t. 

Better  estimates  can  be  made  for  both  the  homogeneous  and 
Inhomogeneous  equations  if  f  is  space -independent .  Existence 
follows  from  the  general  case.   We  write  the  solution 

S|+iLf  =  0  (6. lit) 

as  the  sum  of  its  projection  into  the  nullspace  of  L  and  the 
orthogonal  complem-ent; 

f  -  f  +  f 

f  =  V  p  V   ,    p  =--  {f  ,i'    )  (6.15) 

a 

(f,^J  =   0   . 
The  two  components  satisfy 
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H  +  -  Lf  =  0 


(6.16) 


The  fluid  part,  f.  Is  constant  In  time.   From  (2.13),  we  derive 

||f  II  <  exp(-kLt/e)  ||f^  II  (6.17) 

and  conclude  that  f  decays  exponentially  to  zero  in  a  time  of 
order  e.    Exactly  as  before,  we  extend  this  estimate  to 

nJ  [f]  <  a^  exp(-nt/e)N^  [f^]  ,    r  >  2  (6.l8) 

with  the  same  constant  a  as  previously. 

Next  we  consider  an  inhomogeneous  term  which  is  orthogonal 

to  ii 
a 

|l+-Lf=vg,     f   =0 
dt    e        ^    '  o 

(6.19) 
(vg,  if  J   =   0 

The  same  decomposition  (6. 15)  yields  f  =  0  and  the  a  priori 
estimate 

N^[f]  <   exp(-ix(t-s)/e)N^[  vg(s)  ]ds 

(6.20) 


^ 


o 


<  ^  N^[vg]  . 

This  is  improved  into  a  maximum  estimate  exactly  as  in  the  space- 
dependent  case,  but  we  note  that  e/\i   takes  the  place  of  t   in 
(6.10)  (the  time  t^  is  still  implicit  in  N^); 
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rxS 


Nj  [f]  <  (b^n  +  c^)e  nJ  [g]  ,    r  >  3 


(6.21) 


We  see  that  the  solution  is  smaller  than  the  inhomogeneous  term 
by  a  factor  £  in  this  case 

We  summarize  the  principle  results  in 
Theorem  2:   The  solution  to  the  space -independent  equation 


e 


^  +  T  Lf  =  vg 


\ 


(vg,  ^  )  =  0  ,     (f  ,  ^  )  =  0 


(6.22) 


for  finite  N   [g]  and  n""  [f^].  r  >3,    satisfies 


N^  [f]  <  ^r   exp(-Ht/£)r  [f^] 


+  (h/n  +  c^)£  N^  [g] 


(6.23) 


If  X  occurs  as  a  parameter  in  f  ,  we  merely  replace  N 
by  W^  in  this  theorem. 

Theorem  2  allows  us  to  im.prove  Theorem  1  in  the  space- 
dependent  case  if  vg  and  f  are  orthogonal  to  i^.      We  write 


1    2 
f  =  f  +  f 


(6.2^) 


where 


||i.|Lf^  =  vg,    f^f„ 


(6.25) 


2    12 
Df   +  ^  Lf 


-  e-v 


2 
f   =  0 
o 


(6.26) 


Applying  Theorem  2  to  (6.25),  we  have 


3^  - 


(6.27) 
+  (b/i^  +  c^)  e¥;  [g]  . 

Noting  that  ^^^i  [^-Vf^]  <  W^  [f^].  Theorem  1  applied  to  (6.26) 

yields 

Wj:^  [f2]  <  (b^t^  +  ec^)  ^l   [f^]  (6.28) 

Combining  these  estimates,  we  have 
Theorem  3:   The  solution  to 


Df  +  ^  Lf  -  vg 


(vg,  i/    )    =   0    ,         (r    ,    if    )    =   0 
\  to?   Q,  /         ^  o   a 


(6.29) 


satisfies  the  estimate,  for  r  >  3? 

i^lll    [f]  <  (1  +  b^t^  +  £C^)|a^exp(-iit/e)  ^J  [f^] 

+  (b/^i  +  c^)e¥^  [g]]   .  (6.30) 

If  vg  and  f  are  orthogonal  to  V'  ,  the  solution  rapidly  decays 

to  a  magnitude  of  order  e;  but  to  obtain  this  strong  estimate, 

2  - 
we  lose  a  derivative  of  g  and  a  factor  (l+|  )^. 
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7-   The  Nonlinear  Equation 

We  consider  the  domain  and  boundary  conditions  as  before 

and  Initial  functions  with  a  finite  value  of  W:;  [f  ].   We  look 

p   o 

for  a  solution  of  (2.21)  by  Iteration  on  the  nonlinear  term, 

Df"  +  |Lf"  =  vg^-^  (7.1) 

gn-1  ^  p^^.n-l^^n-1)  (^,2) 

First  we  estimate  #:;  [  p]  from  (5.6).   For  this  purpose  it  Is 

convenient  to  introduce  the  auxiliary  norms 

i_ 
N,  ' 


[V^f]  =||Z  (Ni,[fij))T  (7-3) 


and   sim_llarly  for  ^   .      We   note   the   elementary  Inequalities 
N^iV^f]    <   N^[V®f]    <   3^/^   N^[V^f] 

^^[V^f]   <  i^LV'^f]   <^^^^  »^[V^f]    . 
Taking  norms   in    (5-6)   and  recalling    (4.l6)   and    (4.17), 

<   2/3   T^s^(*^^    [f]f 


(7.4) 


»- 


.[V^T]  <  2Yo^N^ff]  »j^[S/^f]  +  N^[7f]  »^[Vr]j 
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we  obtain 


(7.5) 


y^  =  24  /3  Y^s^    . 
Similarly,  we  derive 

»l  [r(f.g)]  <  Yi  ^l  [f]  ^3  [g]  •  (7.6) 

From  Theorem  1  applied  to  (7-1)  we  have 


^l   [f"]  <a.^^l   [fj  +  t^ll^  [g^-^] 


t-,  -    b^t   +  EC-,     . 

1  30     3 


(7.7) 


The  Iteration  Is  contracting  If  we  choose  J^^  [  f  ]  sufficiently 
small.   Specifically,  if  we  take 

4a^Y^t^  ^ij   [f^]  <  1  ,  (7.8) 

then  all  Iterates  are  bounded 

W^  [f"]  <  2a^  Jf^  [f^]  <  l/(2Y^t^)   .  (7.9) 

Taking  differences  in  (7.I)  and  (7-2) 

^l   [g^--  g"-']  <  Yi^^  [A  f"^-']  w|  [f"-  f"-'] 
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and  we  see  that  the  iterations  contract.   From  the  explicit 
expression  (7-7)  for  t-,  ,  we  conclude  that  (7-8)  will  be 
satisfied  if  t   is  sufficiently  small,  viz., 

t-  <  ^ £T  (7.10) 

where  a  =  (4a^b^Y-,  )    and  t  =  c^b^  depend  only  on  the  col- 
lision operators  L  and  V    and  the  domain  (through  s  ).   A 

bounded  solution  exists  for  0  <  t  <  t   if  t   and  W,    [f    ]    are 

o     o      3   o 

appropriately  related  as  in  (7-10).   Uniqueness  and  continuous 
dependence  on  f  follow  in  the  usual  manner.   We  state 
Theorem  4:   A  unique  solution  to  the  Boltzmann  equation 

Df  +  -  Lf  =  vp(f,f) 

exists  for  a  time  t  which  is  related  to  the  Initial  norm 

o 

W^  [f  ]  by  (7-10).   The  solution  satisfies  the  bound  {1-9), 
¥^  [f]  <  2a-,^^  i^^o-^'  ^^^^'^   ^^'^^  period  of  time. 

Setting  e  =  1  in  this  theorem  yields 
Corollary  1:   A  solution  to 

Df  +  Lf  =  vP  (f  ,f ) 

satisfying  the  bound  ^   [f]  <  2a.^   [f^]  exists  for  a  time 

to  <  -/^l   [^o]  -  ^- 

Similarly,  we  obtain 

Corollary  2:   A  solution  to 

Df  +  Lf  -  evr(f.f) 
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exists  for  t  <  t  /e  where  t   Is  exactly  as  given  In  (7. 10). 
This  result  can  be  obtained  directly  from  the  theorem  by  the 
substitutions  £  — >1,  Y-|  — >  ey-,  . 

A  further  corollary  we  distinguish  as 
Theorem  5:   A  unique  bounded  solution  to 

Df  +  ^  Lf  =  evp(f,f) 

exists  for  t  <  t   =  a/e¥^[f  ]  -  ex   . 

This  Is  obtained  from  Theorem  4  by  the  substitution  y-,  — >  ey-,  . 
All  four  results  can  be  summarized  by  the  statement  that 
a  solution  exists  for  a  sufficiently  small  Initial  perturbation 
from  equilibrium;  and  the  length  of  time  for  which  existence  Is 
proved,  measured  In  collision  times,  Is  approximately  Inversely 
proportional  to  the  relative  size  of  the  Initial  perturbation. 
But  Theorems  4  and  5  are  distinguished  from  the  corollaries  to 
Theorem  4  ty  their  closer  relation  to  tte  linear  equation,  iff  Is  the 
solution  of  the  linear  equation 

Df*  +  I  Lf""  -  0  ,    ^0  =  ^0'  (7.11) 

the  difference 

J25  =  f  -  f*  (7-12) 


satisfies 

where  1  =  0  or  1  under  the  conditions  of  Theorems  4  and  5 
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respectively.   From  Theorem  '^ ,    recalling  that  v|   Is  orthogonal 

to  lA  ,  we  conclude  that 
a 

The  right  side  is  of  order  e  for  the  total  time  of  existence 
referred  to  in  Theorems  4  and  5-   In  Theorem  4,1=0  and 
t   =0(1),  while  in  Theorem  5,  i  =  1  and  t   =  0(l/e).   Thus  f 
differs  from  f  by  an  error  of  order  e  over  the  entire  time  for 
which  the  existence  of  f  is  known.   Under  the  conditions  of 
Corollaries  1  and  2,  the  related  linear  solution  can  only  be 
shown  to  approximate  the  nonlinear  solution  for  a  time  which  is 
short  compared  to  the  proved  time  of  existence;  [o(l)  in 
Corollary  1  and  0(l)  in  Corollary  2].   In  the  one  case  (Theorems 
4  and  5)  the  linearization  is  valid  for  a  time  long  compared  to 
a  mean  collision  time;  in  the  other  case  (corollaries)  only 
for  a  time  comparable  to  a  collision  time. 

The  precise  estimates  of  Theorem  ^5  could  be  used  in  the 
nonlinear  existence  proof  to  improve  Theoremis  4  and  5  (but  not 
the  corollaries).   This  improvement  is  only  by  a  factor  log  l/e; 
i.e.  existence  can  be  shown  for  t^  =  O(log  l/e)  under  Theorem  4 
and  (l/e) log  l/e  under  Theorem  5.   The  improvement  is  so  slight 
because  the  accurate  estim.ate  of  Theorem  3  is  obtained  at  a 
cost  of  reduced  differentiability. 

In  summary,  we  note  that  although  the  convergence  proof 
as  given  is  not  sufficiently  delicate  to  distinguish  those 
cases  (Theorems  4  and  5)  where  the  mean  free  path  of  the  unper- 
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turbed  eq-uilibrlum  co/e  is  small  compared  to  the  length  scale 
of  variation  of  the  perturbation  f,  this  distinction  can  be 
made  by  use  of  Theorem  3- 

We  conclude  this  section  with  a  brief  description  of  the 
existence  theorems  which  can  be  obtained  by  these  methods  for 
the  space -independent  problem  (although  this  is  not  our  pri- 
mary purpose).   First,  since  x  does  not  appear,  we  can  replace 
the   constant  y-,  in  (7-5)  by  y  .   Second  we  may  set  £  =  1 
since  this  only  alters  the  time  scale.   We  obtain  a  solution  of 

II  +  Lf  =  vg 

(7.15) 

g  =  nr,r) 

by  iteration  using  (6.25) 

^^[f]   <  (b^ix  +  c^)N^[g]  (7-16) 

and  (^.l6) 

}^^[r{r^,r^)]   <yJi^[f^]N^[f^]    .  (7-17) 

We  find  that  a  solution  exists  for  all  time  and  staisfies  the 
bound 

N^[f]  <  2a^N^[f^]  (7.18) 

if  N^t^Q]  is  taken  sufficiently  small,  specifically 

2Y^(byn  +  c^)a^N^[f^]  <  1  .  (7.19) 

By  a  slight  modification  of  this  analysis  we  can  show  that 
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these  solutions  are  not  only  bounded  but  decay  exponentially 
to  zero.   First  we  assume  in  the  separation  F  =  co  +  oj^f  that 
o)  is  chosen  with  the  parameters  of  the  ultimate  equilibrium 
(they  are  simply  determined  by  the  total  mass  and  energy  as 
given  by  f^).   Then  (f ,iJ    )  =  0  and  the  approach  to  equilibrium 
consists  in  showing  f  — >  0.   If  we  assume  that  N[e°'^vg]  is 
bounded  uniformly  in  t ,  0  <  a  <  \i,    then  (6.20)  can  be  improved 
to  read  (:Aiith  e  =  l) 

^^'f^^^^f'^  <  -^^   [e^Sg]  (7.20) 

whence  (6.23)  becomes 

N^[e"^f]  <  a^e>:p(-(^i-a)t)N^[f^] 

(7.21) 
+  (by(^i-a.)  +  c^)N^[g]  . 

Exactly  the  same  procedures  as  before  yield  the  result  that  the 
solution  to  (7.15)  satisfies  the  bound 

K%[f]  <  2a^N^[f^lexp(-at)  (7-22) 

provided  that 


Prom  (7.23)  we  see  that  the  decay  exponent  a  can  be  taken 
arbitrarily  close  to  i-l  if  the  initial  value  is  small.   This  is 
not  surprising  since  the  equation  is^  in  this  case^  approximately 
linear  and  |j.  is  the  distance  of  the  spectrum  from  the  origin. 
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For  any  finite  initial  value,  the  ultimate  decay  exponent  for 
large  time  is  arbitrarily  close  to  |j,  since  N^[f]  becomes  small 
ultimately  (7.22),  and  the  theorem  can  then  be  reapplied  with 
a  new  choice  of  a  as  in  (7.25). 

These  spatially  homogeneous  results  can  be  compared  to 
others  of  Carleman  [1]  (also  [6],  Part  l).  Wild  [2]  and 
Morgenstern  [J>] ,   and  Povzner  [^]-   Carleman 's  result  is  for 
hard  spheres  and  an  isotropic  as  well  as  homogeneous  f;  Wild 
and  Morgenstern  consider  only  the  cutoff  Maxwellian  (v  = 
constant);     Povzner 's  equation  is  comparable  to  ours  in 
generality.   Carleman 's  class  of  initial  functions  f   is  some- 
what larger;  Morgenstern ' s  and  Povzner 's  initial  functions  are 
much  more  general.   On  the  other  hand,  Povzner  and  Morgenstern 
prove  nothing  about  boundedness  or  the  approach  to  equilibrium, 
and  Carleman  proves  a  relatively  weak  form  of  the  approach  to 
equilibrium. 


8.   Comparison  with  Fluid  Dynamics 

That  the  solution  of  the  nonlinear  Boltzmann  equation  is 
approximated  by  a  fluid  flow  when  the  mean  free  path  is  small 
is  a  consequence  of  known  results  for  the  linear  equation  (Ref.  [12]) 
using  the  close  correspondence  between  the  linear  and  nonlinear 
equations  that  has  just  been  shown.   For  completeness  of  this 
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presentation  we  shall  give  a  brief  discussion  of  this  theory. 

Under  the  conditions  of  Theorem  h   we  can  show  that  the 
Boltzmann  solution  Is  approximated  within  0(£)  by  a  solution 
of  the  linear  Invlscld  Euler  equations.   The  viscosity  coef- 
ficient Is  of  order  e;  th  s  viscosity  can  be  Ignored  to  this 
order  for  a  finite  time.   Under  the  conditions  of  Theorem  5^ 
the  Boltzmann  solution  Is  approximated  within  0(e)  by  a 
solution  of  the  linear  Navler-Stokes  equation  for  the  whole 
time  for  which  existence  Is  shown.   Viscosity  cannot  be 
ignored  for  such  times  of  order  l/e. 

The  Hllbert  theory,  in  Its  linear  version,  expands  the 
solution  of  (2.22)  as  a  power  series 

fj.  ==  f°  +  ef^  +  ...  ,  (8.1) 

despite  the  fact  that  e  =  0  is  singular.   The  successive  terms 
are  determined  by  an  algorithm  which  Involves  solution  of  an 
integral  equation  for  the  2, -dependence  and  of  a  linear 
fluid-dynamic  Euler  equation  (for  f°)  and  inhomogeneous  Euler 
equations  (later  terms)  for  the  (x,t ) -dependence .   The  Hllbert 
series  (8.1)  is  uniquely  determined  term  by  term,  by  the  initial 
fluid  state  (also  assumed  to  be  a  formal  power  series).  In  other 
words,  to  a  given  fluid  state  the  Hllbert  algorithm  assigns 
a  unique  distribution  function  f^^.   The  truncated  series  (8.1) 
satisfies  an  inhomogeneous  Boltzmann  equation  with  an  inhomo- 
geneous term  that  has  been  estimated  and  shown  to  be  0(e  )  in 
an  appropriate  norm  for  any  finite  time  0  <  t  <  t^.   The  series 
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Is  therefore  asymptotic  to  a  true  solution  of  the  Boltzmann 
equation  which  has  the  same  initial  values. 

The  initial  values  that  can  be  met  by  this  procedure  are 
very  special.   To  generalize  the  Hllbert  theory  a  different 
formal  expansion 

f  (^,x,et)  =  f°  +  £f^  +  ...  (8.2) 

can  be  found  which  also  satisfies  the  Boltzmann  equation  term  by 

term.   But,  whereas  f„  is  uniquely  determined  by  the  fluid  state, 

n 

f   Is  uniquely  determined  by  its  orthogonal  complement,  and  each 
M- 

term  in  the  expansion  of  f  decays  uniformly  as  exp( -M-t/e) . 

Although  f„  and  f   are  determined  by  the  projections  of  f   into 

the  nullspace  of  L  and  its  complement,  the  initial  f^   and  f,, 

are  not  just  the  projections  of  f  .   By  summing,  f  =  f.,  +  f  , 

O  n     M- 

essentially  arbitrary  initial  data  can  be  matched,  and  we  can 
find  an  appropriate  f„  (satisfying  a  suitably  modified  initial 

n 

fluid  state)  which  is  asymptotic,  after  a  short  transient,  to 
general  smooth  solutions  of  the  linear  Boltzmann  equation. 

The  Chapman -Enskog  procedure  is  an  alternative  to  Hllbert, 
one  which  uses  the  same  terms  (8.1)  as  functions  of  ^  but 
Inserts  a  modified  (x,t)  dependence.   Instead  of  Euler  and 
inhomogeneous  Euler  equations,  the  fluid  parameters  in  f„  are 
determined  by  Euler,  Navler-Stokes,  and  successively  highly 
order  differential  equations. 

This  expansion  is  also  asymptotic  to  a  solution  of  the 
Boltzmann  equation  provided  that  one  can  suitably  estimate  the 
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solutions  to  the  fluid  equations.   The  Hllbert  theory  has  the 
advantage  that  It  can  be  carried  to  arbitrarily  high  order  In  e 
by  solving  hyperbolic  differential  equations  which  have  a  known 
theory.   For  the  Chapman-Enskog  procedure  there  Is  no  theory 
beyond  the  second  step  (Navler-Stokes) .   But  the  (linear) 
Navler-Stokes  solution  can  be  bounded  uniformly  for  all  time, 
whereas  the  equivalent  Kllbert  solution  grows  linearly  with 
time.   Thus  the  second  Chapman-Enskog  term  Is  asymptotic  uni- 
formly In  time;  for  our  application  we  need  Its  validity  for 
0  <  t  <  tjz. 

We  shall  sketch  the  proofs  that  the  first  order  Hllbert 
approximation  Is  asymptotic  to  the  linear  solution  under 
conditions  appropriate  to  Theorem  4,  t  =  0(l),  and  the  second 
order  Chapmian-Enskog  solution  Is  asymptotic  to  the  linear 
solution  under  conditions  appropriate  to  Theorem  5?  t  =  0(l/£); 
for  more  complete  results  see  Ref.  [12]. 

We  take  an  initial  function  f   of  suitable  regularity  and 
construct  the  solution  f  of  the  linear  equation 

Df  +  I  Lf  =  0   .  (8.3) 

For  the  first  case,  covered  by  Theorem  4,  we  construct  a 
locally  Maxwelllan  comparison  function 

a 
The  p   are  related  to  the  (linearized,  dlmensionless)  nominal 

fluid  state  by 
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f      Po   =   P 

<    Pi  =  u. 


(1  =  1,2,3) 


(8.5) 


p^  =  (V2)2t 


We  specify  the  p   In  f  as  functions  of  x  and  t  by  choosing 
them  to  be  solutions  of  the  (linearized,  dlmenslonless)  Euler 
equations 


1^  +  dlv  u  =  0 

^"l  ^  ap   ^  BT 
dt  dx.         dx. 


3  St 

2  ^ 


+  dlv  u  =  0 


(8.6) 


with  Initial  data  taken  from  f 


o 


p  (x,0)  =    {ii    ,f    ) 
^a^         '         ^  a   o' 


(8.7) 


The  difference 


^   =   f    -   f 


E 


(8.8) 


satisfies   the   equation 


D/5  +  -  L/5  -    -  Df 


E 


(8.9) 


By  the  choice  of  Initial  data  we  have 


(^■*a: 


0 


(8.10) 


and,  from  the  fact  that  p   In  f_  satisfies  the  Euler  equations, 

Ot      E 
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we  can  easily  verify  that 

(Dfg,  ^^)  =  0  .  (8.11) 

Inspection  of  Theorem  3  shows  that,  after  a  short  exponential 
transient,  f)   becomes  small  of  order  e.   We  note  that  It  Is  the 
same  estimate  (Theorem  3)  which  shows  that  the  linear  solution 
approximates  the  Euler  solution  f_  and  also  that  the  nonlinear 
solution  approximates  the  linear  solution. 

For  the  second  case,  covered  by  Theorem  5,  we  choose  a 
slightly  m.ore  complicated  comparison  function 

As  before,  f  Is  locally  Maxwelllan 

f AT  =  51  P  ^  (8-13) 

a 

but  the  p   are  now  taken  to  satisfy  the  linear  Navler-Stokes 

equations 

/  1^  +  dlv  u  =  0 

^  +  |£_  +  1^  =  eTi[Au.  +  i  ^  (dlv  u)]       (8.14) 
dt     dx .    dx .     '    1    5  ox .  ^      '  ^  ^  ' 

^  -^  +  dlv  u  =  e  A  AT 
2  dt 

(sT]  and  eA  are  the  dlmenslonless  viscosity  and  heat  conductivity 
coefficients).   By  inspection  of  (8.12)  one  should  expect  f^  to 
approximate  f  to  within  0(e  ).   This  Is  Indeed  possible  for  a 
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time  which  Is  0(l).   To  accomplish  this  the  p   In  f  would 
have  to  be  assigned  appropriate  Initial  values  which  are 
obtained  from  f  by  a  rather  complex  computation  (see  Ref. 
[12]).   But,  since  we  are  Interested  In  times  of  order  l/s, 
we  shall  only  be  able  to  show  that  f  approximates  f  to  within 
0(e).   This,  as  it  turns  out,  will  permit  us  to  use  the 
simpler  initial  values  (8.7).   Summarizing,  we  specify  the 
leading  term  f,,  by  (8.13)  where  the  p   are  solutions  of  the 
Navler-Stokes  equations  with  given  initial  data  (8.7). 

The  second  term  in  the  comparison  function  f,,  is  obtained 
by  the  Chapman-Enskog  algorithm.   Formal  substitution  of 

P  =  f  -  f jj  (8.15) 

into  the  Boltzmann  equation  yields 

■Djd  +  J  Lj6  =    -   Dfj^  -  Lfj^  (8.16) 

We  would  like  to  specify  f„  to  satisfy 

Lf^  +  Df^  =  0  (8.17) 

since  this  would  make  the  inhom.ogeneous  term  in  (8.16)  of  order 
e.   This  is  an  Integral  equation  for  f  ,  with  Df.  a  known 
inhomogeneous  term.   But  the  nontrlvlal  nullspace  of  L  requires 
a  compatibility  condition  on  Df,,  which  is  not  satisfied.   From 
our  earlier  result  (8.11)  we  know  that  compatibility  would  be 
satisfied  if  p   satisfied  the  Euler  equations.   Inspecting 

a 
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and  replacing  bp  /ht   by  their  Navier-Stokes  equivalents  fr 


a 


om 


(8.l4)  we  can  set 

Dfjj  =  5^  +  e6^  (8.19) 

where  In  5  we  collect  the  Euler  terms  from  (8.l4)  and  in  5 
the  viscous  and  heat  conducting  terms.   It  is  easily  verified 
exactly  as  in  (8.11),  that 

(^^,  5J  =  0  (8.20) 

We  therefore  choose  f,-  not  as  the  (nonexistent)  solution  of 
(8.17),  but  as  the  unique  solution  of 

Lfj^  +  60  =  0  (8.21) 

This  leaves  a  formal  remainder  in  (8.16)  of  order  e, 

DJ2^  +  I  L/^  -  -  eDf^  -  65^  .  (8.22) 

We  note  that  the  separation  (8.19)  is  strictly  formal;  the  p 
as  solutions  of  the  Navier-Stokes  equations  are  complicated 
functions  of  e  and,  in  general,  are  non-expandable.   But  the 
remainder  in  (8.22)  can  be  estimated  to  be  0(e),  given  suffi- 
ciently smooth  initial  data  for  the  Navier-Stokes  solution. 

The  solution  of  (8.21)  for  f^,  is  an  explicit  linear  com- 
bination of  velocity  and  temperature  gradients  with  coefficients 
which  are  functions  of  I  obtained  by  solution  of  the  integral 
equation  (which  can  be  transformed  into  a  Fredholm  equation 
with  compact  kernel).   It  is  easily  verified  that 
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{■^^,    Df^)    =   0   .  (8.23) 

This  is  essentially  the  statement  that  p   satisfy  the  Navler- 
Stokes  equations [cf.  the  missing  orthogonality  condition  for 
(8.17),  (I/'  ,  ^^iq)   /  0];  this  is,  in  fact,  the  kinetic  derivation 
of  the  Navier-Stokes  equations.   ¥e  are  now  in  a  position  to 
refer  again  to  Theorem  3-   Since  the  Inhomogeneous  term  in 
(8.22)  is  of  order  e  (uniformly  for  all  time  by  Navler-Stokes 
estimates)  and  is  orthogonal  to  i'   ;and  since  the  initial  value  j6^ 
is  orthogonal  to  i/'  ,  we  conclude  that  j6   is  0(e)  for  t  =  0(l/e) 
[we  do  not  obtain  0(e  )  because  of  the  factor  t   =  0(l/e)  in 
(6.30)].^ 

The  solutions  of  the  nonlinear  Boltzmann  equation  obtained 
in  the  previous  section  have  been  shown  to  approximate  fluid 
behavior  only  to  the  extent  that  they  are  also  linear.   While 
these  fluid  comparisons  can  be  improved  to  some  extent,  they 
cannot  be  greatly  extended  by  our  present  methods.   For 
example,  the  existence  theory  is  basically  unsuited  to  des- 
cribe a  finite  strength  shock  wave.   A  shock  wave  of  finite 
amplitude  would  have  a  thickness  0(e)  and  would  require  handling 
estimates  N[Vf]  -  0(e~"^/^),  N[V^f]  =  0(e~^/^),  etc.   A 
Sobolev  estimate  would  not  be  strong  enough  to  rule  out  large 


The  stronger  result  quoted  in  Ref . [12] (p. 171 )  is  in  error. 

Although  the  linear  Navier-Stokes  solution  in  a  bounded 

domain  decays  exponentially,  the  decay  is  slow  as  exp(-et) 

00 

Terms  such  as  /  f  dt  will  produce  a  factor  l/e  just  as  it 

o    , 
arises  from  t  ^^  l/e. 
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excursions  In  f  itself. 

Under  Theorem  4,  a  fully  developed  shock  wave  with 

strength  of  order  e  would  have  a  finite  thickness,  0(l). 

2 
Uncer  Theorem  5,  with  nonlinearity  of  order  e  ,  the  fully- 
developed  thickness  would  be  0(l/e).   Thus  weak  shocks,  and 
even  much  steeper  profiles,  are  contained  within  the  theory. 
However,  the  time  required  to  observe  the  development  of  a 
steady  limiting  shock  profile  would  be  0(l/e)  under  Theorem 
4  and  0(l/e  )  under  Theorem  5.   We  could  therefore  observe 
only  the  initial  stages  of  shock  formation. 

One  may  hope  that  stronger  nonlinear  results  for  the 
Boltzmann  equation  will  be  obtained  by  making  efficient  use 
of  nonlinear  Navier-Stokes  estimates. 
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Appendix:   Estimate  of  the  Quadratic  Collision  Integral-^ 

In  a  more  complete  notation,  the  term  Q  in  (2.1)  can  be 
written 

Q(P,F)  -  /(P'p'  -  FF,  )dD. 

^  ^  (A.l) 

dQ   =   B(0,V)d0d£d^^  . 

Introducing   f  as   in    (2.4),   we    obtain  the   bilinear  form 

vr(f.g)    =/i(f'g|   +   fjgi    -   fgi    -   fig)a3jdfi  (A. 2) 

where 

V   =  /cD^dn    .  (A.  3) 

To  estimate    f"""  we  write 

IvFI    <   vFi  +   vPg 

vFi   =  /ilfgi  +  fighjdfi  [kA) 

v\~'^   =  /i|f'g^  +   fjg^|a)2^fi    . 

To  estimate   F-,    ^e   need   the  bound 

o   =  /co|dfi  <  a    V  (A.  5) 

for  some  constant  a    .      This  is  easily  obtained  from  the 
explicit  forms 

v(^)  =  2-n-/BQ(h-e|  )a)(Ti)dri 

,  (A. 6) 

a(^)  =   27r/B^(h-^|)a)2(ri)dii 


9.   For  the  notation,  transformations,  and  supporting  estimates 
used  in  this  Appendix,  we  refer  to  Ref.  [9]- 
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where 


B^(V)  =  /B^(e,V)d9   .  (A. 7) 


The  bound  (A. 5)  follows  from  (A. 6)  as  an  elementary  exercise 
if,  for  example,  B  (V)  has  a  uniformly  bounded  derivative 
(this  is  true  for  all  hard  power  law  potentials). 
From 

|f|  <  (l+^2)-^/2  j^^ffj  ^  jj^[^] 


we  obtain 


<  (1+^^)"^/^  N  [g]  <  N  [g] 


ilfgl  +  f^gl  <  {l+i^)~''^^   N^[f]N^[g] 


and  conclude  that 


^Fi  <  (l+e^)-^/^  a^vN^[f]N^[g] 


or 


N^Lrj  <  «oN^[f]N^[g]  .  (A. 8) 

To  estimate  Vo'   ^®  transform  the  integral  as  in  Ref.[9] 
(here  we  use  R(v,w)  for  the  kernel  Q(v,w)  in  Ref.  [9]) 

vfs  =/i|f(^+v)g(^+w)  +f(^+w)g(^+v)|u32(^+v+w)5-^^^  dv  dw. 

V 

(A. 9) 
We  recall  that  w  is  integrated  over  the  plane  perpendicular  to 
V,  and  then  v  is  integrated  over  the  full  3-dimensional  space. 
Resolving  ^  into  components  4.  and  |„  parallel  to  v  and  w 
respectively,  we  have 
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co^(4+v+w)  =  CD2(^-^+v)a3^(?,2+w)  (A. 10) 

From  eq.  (6o)  of  Ref.  [9]  we  have 

J   co2(f,2+w)R(v,w)dw  <  Q,  V  .  (A. 11) 


From 


[l+{i+vf][l+{i+^f]   >   l+(^+v)^  +    (^+w)2 


1+2,     +     (^+v+w) 


we    obtain 


>   1+^^ 


f(l+v)g(^+w)|    <    (l+^^)"''/^N^[f]N^[g]     .  (A. 12) 


Substituting  (A. 11)  and  (A. 12)  Into  (A. 9),  we  find 

vPg  <  Q^(l+l^)"'^/\[f]N^[g]  /^CD*(^^+v)dv  .       (A. 13) 

The  last  Integral  Is  bounded  (it  attains  its  maximum ;,  4(27r)  ^     , 
at  I-,  =  O),  thus  we  have 

Po  ==  MStt)^/^  Q^  . 
Finally  we  have  the  desired  result, 

where  v^  =  a  + 6  •   A  trivial  modification  of  the  above  argu- 
'0    o   o 

ments  yields  the  companion  estimate 
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